Influence of polymer-pore interactions on translocation 
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We investigate the influence of polymer-pore interactions on the translocation dynamics using 
2D Langevin dynamics simulations. An attractive interaction can greatly improve translocation 
probability. At the same time, it also increases translocation time slowly for weak attraction while 
exponential dependence is observed for strong attraction. For fixed driving force and chain length 
the histogram of translocation time has a transition from Gaussian distribution to long-tailed distri- 
bution with increasing attraction. Under a weak driving force and a strong attractive force, both the 
translocation time and the residence time in the pore show a non-monotonic behavior as a function 
of the chain length. Our simulations results are in good agreement with recent experimental data. 

PACS numbers: 87.15.Aa, 87.15.He 



The transport of a polymer through a nanoporc plays 
a critical role in numerous biological processes, such as 
DNA and RNA translocation across nuclear pores, pro- 
tein transport through membrane channels, and virus in- 
jection. For a polymer threading through a nanopore, 
loss of available configurations due to the geometric con- 
striction leads to an effective entropic barrier. Kasianow- 
icz et al. [l[ demonstrated that an electric field can drive 
single-stranded DNA and RNA molecules through the 
water-filled a-hemolysin channel and that the passage of 
each molecule is signaled by a blockade in the channel 
current. These observations can directly be used to char- 
acterize the polymer length. Due to various potential 
technological applications [3, 0] , such as rapid DNA se- 
quencing, gene therapy and controlled drug delivery, the 
polymer translocation has become a subject of intensive 
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As to translocation, one of the basic questions concerns 
the dependence of the translocation time r on the system 
parameters such as the polymer chain l eng th N 0, @, Hoi 
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Dore length L and 
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and secondary structure [3j, |4j, [6|, [22 
pore width W [ljj[ , driving force F 
23l . |24j , and polymer-pore interaction 

In a recent experiment, 0,0 striking differences were 
found for the translocation time distribution of poly- 
deoxyadenylic acid (poly(dA)ioo) and polydcoxycytidylic 
acid (poly(dC)ioo) DNA molecules. The origin of the 
different behavior was attributed to stronger attractive 
interaction of polyfdA) with the pore. Also, recently 
Krasilnikov et al. [9( have investigated the dynamics of 
single poly (ethylene glycol) (PEG) molecules in the a- 
hemolysin channel in the limit of a strong attractive 
polymer-pore attraction. The result for the residence 
time in the channel shows a novel non-monotonic behav- with a cutoff of 2.5a and interaction strength e pm . This 



ior as a function of the molecular weight. 

On the theoretical front, not only the quantitative 
but also the qualitative picture of the polymer-nanoporc 
interactions is still elusive. Based on a Smoluchowski 
equation with a phcnomenological microscopic poten- 
tial to describe the polymer-pore interactions, Lubcn- 
sky and Nelson [3] captured the main ingredients of the 
translocation process. However, when comparing with 
experiments, their model is not sufficient. Numerically, 
Tian and Smith (23[ found that attraction facilitates 
the translocation process by shortening the translocation 
time, which contradicts experimental findings [1, 

To this end, in this letter we use Langevin dynam- 
ics (LD) to investigate the influence of polymer-pore in- 
teractions on translocation. In our numerical simula- 
tions, the polymer chains are modeled as bead-spring 
chains of Lcnnard- Jones (LJ) particles with the Finite 
Extension Nonlinear Elastic (FENE) potential. Ex- 
cluded volume interaction between monomers is mod- 
eled by a short range repulsive LJ potential: ULj(r) = 
- (f) 6 } + e for r < 2 1 / 6 a and for r > 2 1 / 6 a. 
Here, a is the diameter of a monomer, and e is the 
depth of the potential. The connectivity between neigh- 
boring monomers is modeled as a FENE spring with 
Ufene{t) = — ^At-Rq ln(l - r 2 /i?p), where r is the dis- 
tance between consecutive monomers, k is the spring con- 
stant and Rq is the maximum allowed separation between 
connected monomers. 

We consider a 2D geometry as shown in Fig. [TJ where 
the wall in the y direction is described as stationary par- 
ticles within a distance a from each other. The pore 
of length L and width W in the center of the wall 
is composed of stationary black particles. Between all 
monomer-wall particle pairs, there exist the same short 
range repulsive LJ interaction as described above. The 
pore-monomer interaction is modeled by a LJ potential 



2 




-14 -12 -10 



FIG. 1: A schematic representation of the system. The pore 
length L = 5 and the pore width W = 3 ( See text for units) 



interaction can be either attractive or repulsive depend- 
ing on the position of the monomer from the pore par- 
ticles. We have also performed numerical calculation 
for the case of a pure short range repulsive LJ poten- 
tial for the pore-monomer interaction. As expected, the 
results for the long range LJ pore-monomer interaction 
approaches that for the pure repulsive pore-monomer in- 
teraction in the limit e pm — > 0. In the Langevin dynam- 
ics simulation, each monomer is subjected to conserva- 
tive, frictional, and random forces, respectively, with [2H 
mvi = -V(Ulj + Ufene) + F cxt - + Ff , where m 
is the monomer's mass, £ is the friction coefficient, 
is the monomer's velocity, and is the random force 
which satisfies the fluctuation-dissipation theorem. The 
external force is expressed as F oxt = Fx, where F is the 
external force strength exerted on the monomers in the 
pore, and x is a unit vector in the direction along the 
pore axis. 

In the present work, we use the LJ parameters e and 
a and the monomer mass m to fix the energy, length 
and mass scales respectively. Time scale is then given 
by thj = (mcr 2 /e) 1 / 2 The dimensionless parameters in 
our simulations are Rq = 2, k/m = 7, fc^T = 1.2, and 
£/m = 0.7. For the pore, we set L = 5 unless otherwise 
stated. A choice of W = 3 ensures that the polymer 
encounters an attractive force inside the pore. We have 
checked that a choice of W — 4 yields similar results. 
The driving force F is set between 0.5 and 2.0, which 
correspond to the range of voltages used in the experi- 
ments [H, [H . The Langevin equation is integrated in time 
by a method described by Ermak and Buckholtz 2(i\ in 
2D. Initially, the first monomer of the chain is placed in 
the entrance of the pore, while the remaining monomers 
are under thermal collisions described by the Langevin 
thermostat to obtain an equilibrium configuration. The 
translocation time is defined as the time interval between 
the entrance of the first segment into the pore and the 
exit of the last segment. Typically, we average our data 
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FIG. 2: (a) The translocation probability as a function of the 
attractive strength for different driving forces, (b) The distri- 
bution of translocation time for different attractive strengths 
under the driving force F = 0.5. The chain length N = 128. 
The data point at e pm = corresponds to a pure repulsive 
pore-monomer interaction 



over 2000 independent runs. 

The translocation probability, Pt rans , is calculated 
as the fraction of runs leading to successful transloca- 
tion at given conditions. Fig. 2(a) shows Ptrans as a 
function of e pm for = 128 under different driving 
forces. Specifically, the numerical results clearly show 



two different regimes. With increasing e pm , Ptrans in- 
creases rapidly first, and then slowly approaches satura- 
tion at larger e pm . It is known that attractive interac- 
tion with the channel can facilitate the translocation of 
metabolite molecule across cellular and organelle mem- 
branes (2^ |. Here, our results show that polymer translo- 
cation through the attractive nanopore shares the same 
character. 

Experimentally, Mcllcr et al. 0, @1 have investigated 
the translocation of homepolynuclcotides of different 
bases: poly(dA)ioo and poly(dC)ioo- The translocation 
time distributions in both cases are well approximated by 
fast-growing Gaussian for translocation time lower than 
the most probable value t p and falling exponentials for 
translocation time larger than t p . The decay time scale 
for poly(dA) is found to be much longer than poly(dC), 
by a factor of ~ 7. There exists also a large difference 
between the value of t p for the two, corresponding to 
1.2 /xs/base for poly(dC), and 3.3 [is/base for poly(dA). 
These differences have been attributed to the base spe- 
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FIG. 3: Translocation time as a function of the attractive 
strength for different driving forces. The chain length N — 
128. Here, e pm = corresponds to a pure repulsive pore- 
monomer interaction 
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FIG. 4: Translocation time as a function of the chain length 
for e P m = and e pm = 1 with F — 0.5. The insert shows 
results for e vm = 3 with F = 1.0 and F = 2.0. 



cific nucleotide-pore interactions, with the adenines hav- 
ing a stronger attractive interaction with the pore as com- 
pared with cytosines. 

In our numerical results, we have found that indeed for 
e pm =2 and 3, the attractive potential has a marked im- 
pact on the shape of the histogram of the translocation 
time as shown in Fig. 2(b). The shape of the histogram 
changes from a nearly Gaussian below the most proba- 
ble value to a long exponential tail. The value of t p as 
well as the characteristic decaying time scale increases 
with e pm . These findings are in excellent agreement with 
the experimental observation of Mcller et al. 0, Q , and 
provide further support that the base specific interaction 
with the pore plays a pivotal role in the translocation 
dynamics of single- stranded DNA and RNA molecules. 

Fig. 3 shows the calculated tF as a function of e pm 
for A — 128 under different driving forces. Initially, 
r increases very slowly with increasing e pm . Then, it 
crosses over to a different regime and increases sharply 
to the asymptotic behavior tF ~ e L / k BT _ rpj^ crossover 
threshold value of e pm increases with increasing F. Sur- 
prisingly, previous numerical work [23j failed to capture 
the essential feature that the translocation time increases 
with increasing attractive base-pore interaction. 

An important element of our analysis is the fact that 
the translocation time can be written as r ~ t\ + T2 + T3, 
where t\, t-i and T3 correspond to initial filling of the 
pore, transfer of the DNA from the cis side to the trans 
side, and finally the emptying of the pore, respectively. 
In the presence of the attractive pore-monomer interac- 
tion and driving force across the pore, n << T2, T3, 
while T2 increases monotonically with TV. For strong 
attraction and intermediate values of A, r is deter- 
mined mainly by T3 related to the emptying of the 
pore. This process involves a free energy difference of 
AF = L(e pm - Fa/2 - /(A)) between the final and the 
initial state. The term /(A) here accounts for the en- 



tropic driving force which should kick in at larger val- 
ues of A and eventually saturate for very long polymers. 
For the region of weak attraction below the threshold, 
AF < and the translocation time depends weakly on 
e pm . Above the threshold when AF > 0, the process 
is activated with a barrier ~ AF and increases rapidly 
with increasing strength of attraction e pm . This accounts 
for the observed crossover behavior of r as a function of 
e pm . In the weak attraction non-activated region, the 
overall r is determined mainly by T2 and its dependence 
on the driving force scales as -F -1 which comes from the 
velocity dependence on F. However, once one enters the 
activated region, the force F also influences the activa- 
tion barrier besides affecting the prefactor and t drops 
off with increasing F much faster than the simple F _1 
behavior. 

For F = 0.5 and a pure repulsive pore-monomer inter- 
action, we have shown in our earlier work [H H3 that 
t ~ N 2lJ for relatively short chains and crosses over to 
t ~ N 1+u for longer chains as shown in Fig. 4, where 
the Flory exponent v = 0.75 in 2D [28|, and the crossover 
length A c ~ 200. The scaling behavior for attractive in- 
teraction strength e pm = 1 is very similar to the pure 
repulsive case. For e pm = 2, we found that A c ~ 310. 
For stronger attractive strength e pm = 3, only r ~ A 2l/ 
is observed for the A values studied under F = 1 and 
F = 2, with no indication of crossover behavior as shown 
in the insert of Fig. 4. 
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weak driving force F = 0.5, the translocation time r has 
a qualitatively different dependence on A as compared 
with the pure repulsive or weak attractive pore interac- 
tion. It has a novel non-monotonic behavior with a rapid 
increase to a maximum at A ~ 14, followed by a decrease 
for 14 < A < 32and an increases again for A > 32 as 
shown in Fig. 5(a). This can be understood by consid- 
ering the different N dependence of n, t-i and T3 in the 
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FIG. 5: (a) Translocation time r as a function of the chain 
length for e pm = 3 and F — 0.5. (b) Residence time r r as a 
function of the chain length for e pm = 3 and F = 0. 

strong attraction limit. For small and intermediate val- 
ues of N, t is dominated by T3. Here the entropic factor 
f(N) in the barrier for T3 fights against the simple power 
law increase in the prefactor accounting for the number 
of monomers needed to cross the pore. This leads to 
an initial increase of r to a maximum value followed by a 
subsequent decrease. Eventually, for larger N, the T2 pro- 
cess (which approaches N 2u asymptotically) takes over, 
leading to the increase of r with increasing A*" again. 

For this case, we found that there is about 20% of the 
total translocation processes in which the polymer enters 
and reexits the cis side of the pore. It is useful to de- 
fine an additional residence time the weighted sum 
of the translocation time and the return time which cor- 
responds to the experimentally measured blockage time. 
For the case with no external driving force, the translo- 
cation probability is very small and the residence time 
is almost all due to return events. Wc have calculated 
the residence time r r for F = and s pm = 3 and the 
result is shown in Fig. 5(b). The N dependence here is 
again non-monotonic similar to the translocation time for 
F = 0.5 except for the absence of the eventual increase 
at the large N limit, due to the absence of the T2 contri- 
bution for the return process. Our numerical result of r r 
for F = is in good agreement with recent experimen- 
tal data of Krasilnikov et al. 0] in which the residence 
time of neutral PEG molecule in a-Hemolysin pore was 
measured. 

To summarize, wc have investigated the influence of at- 
tractive polymer-pore interactions on the translocation 
dynamics via numerical simulation studies of a simple 
course grained model. Our results are in good agreement 
with recent experimental data for driven translocation of 
poly(dA)and poly(dC) molecules, and for the blockage 
time study of poly(ethylene glycol) molecule through a- 
Hemolysin pore. They clearly demonstrate the important 
role of polymer-pore interaction factor in the transloca- 
tion dynamics. 
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